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Abstract. Consider a second order hyperbolic linear system with time inde- 
pendent coefficients whose solution and source term for each time moment are 
sections of a real smooth vector bundle over a closed connected smooth man- 
ifold. Assume that the energy of the system is conserved. Moreover, assume 
that the manifold together with its topological and differential structures and 
a smooth measure on it are known and we are given sets of sources which are 
known to be localized in space and time and dense in some function spaces. By 
measuring how much energy is needed to produce different combinations (i.e. 
sums) of these sources, we reconstruct a vector bundle up to an isomorphism, 
the Riemannian structure on it and the system. We also show that these sets 
of sources are generic in the sense that they can be almost surely constructed 
by taking some sequences of realizations of suitable independent identically 
distributed Gaussian random variables. 



1. Introduction 



In this paper we consider a smooth real vector bundle V with a Riemannian 
structure over a compact closed connected smooth manifold M and a second 
order hyperbolic initial value problem 

+ A{x, d^))u{x, t) = F{x, t), X e M, t > -r, 
u{x, t) |t=_r = 0, dtu{x, t) = 0, 

with some large r > and supp {F) G M x (— r, -(-oo). Here for each time 
moment t G M the solution u{- ,t) = u^{- ,t) : M ^ V and the source term 
F{- ,t) : M ^ V are sections of the vector bundle V. 

Let an elliptic operator A = A{x,dx) be such that the energy of the wave 
produced by the source F, defined by formula 

EAiF, t) = l [ {{dtu'^ix, t), dtu^ix, t)), + {Au'^ix, t), t))x)dV{x), 
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is conserved, that is, if F = for all t>to then Ea{F, t) — Ea{F, to). Here (•, ■)x 
is the inner product in the fiber Tr~^{x) of the vector bundle V assigned by the 
Riemannian structure and dV{x) is a smooth measure on M. 

Assume that the manifold M together with its topological and differential struc- 
tures and a smooth measure dV{x) on it are known. The inverse problem studied 
in this paper is to reconstruct the vector bundle V, the Riemannian structure 
on it and the hyperbolic system by measuring energy of waves u^^ , as functions 
1 1— > E^iFj, t), for certain family of sources Fj, j — 1,2,?,, 

Since the vector bundle is not known, we encounter immediately the problem how 
to specify the sources Fj which we use in the measurements. To approach this 
problem, we assume that in some method we can produce the families of sources 
of the following type: 

(i) C°°(M X /) -dense sequence of sources Fj, j — 1,2, on M x I where 
/ C M is an open interval; 

(ii) C{U X /)-dense sequence of sources Fj, j — 1,2, . . . , on U x I where U is 
an open set on M. 

These families of sources are generic in the sense that they can be almost surely 
generated by taking some sequences of realizations of suitable independent iden- 
tically distributed Gaussian random variables, see Section 3. Roughly speaking, 

we assume that we can just produce a family of physical sources supported in 
the set U X I but we do not know what the functions describing these sources 
are. We call such produced sources the basic experiments. We also assume that 
we can combine disjointly supported basic experiments and delay them in time 
and measure how much energy is needed to perform the combined experiments. 
Using this data, we want to determine the physical model, that is, the bundle 
and the coefficients of the hyperbolic system. 

Similar problems of recovering the operator A{x, d^) and the Riemannian man- 
ifold from energy measurements have been encountered in the inverse boundary 
value problems. The prototype inverse problem is the inverse conductivity prob- 
lem, called also the Calderon's inverse problem (see [6]). Consider a conductivity 
equation 

V-a{x)Vu{x) = inD, (1.1) 
u\dD = /, 

where D C is smooth domain, (j{x) is a positive-definite matrix corresponding 
to the conductivity at the point x, u is the electric potential having boundary 
value /. The inverse conductivity problem means the determination of cr from 
measurements done at the boundary dD. One possibility to define the boundary 
measurements is the following (there are many equivalent formulations): assume 
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that for all boundary values / we measure 



This is equivalent to measuring the power needed to keep the voltage u at the 
boundary dD to be equal to f{x). Physically, the boundary value u\q£, = f 
is produced by electrodes attached on the boundary of the body, and Qa{f) 
corresponds to the power of the heat produced by the caused currents in the 
domain D. Because of concervation of energy, the power produced by the heat is 
equal to the power needed to keep the electrodes at the specified voltages. When 
a is equal to a scalar function j{x) times the identity matrix, the conductivity 
is called isotropic. In this case the measurements are shown to determine the 
conductivity - the problem was solved for smooth conductivities in [24, 20] and 
in the form posed by Calderon in [6] for L°° conductivity in two dimensions in 
[2], see also [5, 13, 16]. 

When the conductivity a in (1.1) is matrix valued (physically, the conductivity 
is anisotropic), it has been been shown that the conductivity a can not be deter- 
mined uniquely and that the best one can hope is to determine cr up to a change 
of coordinates. This has been shown in dimension 2 in [3, 22, 23]. In dimensions 
n > 3 determining the conductivity is equivalent to determining the isometry 
type of a Riemannian metric corresponding to the conductivity, and this is a 
longstanding open problem. For related works, see [10, 11, 14, 18, 19]. 

As another example of inverse problems on a Riemannian manifold which moti- 
vates our study, let us consider the wave equation. 



on a compact Riemannian manifold (M, g) with boundary dM. Assume that we 
know the boundary dM and can measure for all / G C^{dM x R_,_) the final 
energy Qg{f) of the waves u produced by /, that is. 



By energy concervation, this can be considered as the total energy needed to 
force the boundary value u\qmxM.+ to be equal /. When Qg{f) is known for all 
/ G C^{dM X M+), then the isometry type of the Riemannian manifold {M,g) 
can be recoved [15], see also [1, 14]. 

In modern physics, phenomena are often modeled by PDE systems on sections of 
vector bundles. The vector bundles are encountered from the early development 
of the relativistic quantum mechanics to modern quantum field theory. Hence 
the study of inverse problems on general vector bundles brings the mathematical 





(1.2) 
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study of inverse problems closer to modern physics. A result in this direction is 
[17] where the inverse problem for the Dirac equation on a general vector bundle 
over a compact Riemannian manifold with nonempty boundary is studied. 

Let us consider also other examples, analogous to the problem studied in this 
paper, where reconstruction of the bundle structure is interesting. Consider the 
collision of elementary particles. When such collisions are produced, the precise 
parameters determining the collision are not usually known, but one just observes 
that a collision has happened and its consequences can be also observed. Such 
non-exactly known experiments are analogous to the measurements considered in 
this paper. By observing behavior of electrons and other particles the physicists 
deduced the existence of the spin structures, which implies that vector bundles 
with spin structures are important in the large scale models in a curved universe. 
A priori, before the physical observations, it is not clear what the bundle struc- 
ture is. There arc physical models for yet unobserved particles, such as Dirac's 
magnetic monopole, which are based on the use of non-trivial bundles. Thus, if in 
the future one would observe a new particle behaving like a magnetic monopole, 
and the results of the physical experiments for that would coincide with those 
predicted by Dirac's magnetic monopole model, one should ask if there are other 
models giving the same physical predictions. The reconstruction problem of the 
bundle structure studied in this paper can be considered as a highly idealized 
model problem for such a physical problem. 

Reconstruction of unknown "bundle structure" appears also in much more mun- 
dane problems than in the high energy physics. When making observations on a 
chemical system where reactions happen in several compartments (e.g. in differ- 
ent cells and blood in a metabolic model [7, 8]), we can often observe only part 
of the chemical concentrations and only at some compartments. The diffusion- 
reaction equation can be considered as an equation on a vector bundle, where 
the structure of the bundle corresponds to the knowledge, how many substances 
are present in the reactions, and the differential equation corresponds to the 
knowledge what kind of reactions happen in each compartment. Similarly, in 
mathematical finance, in modeling prices on a stock market, one can look prices 
as observations and ask how many hidden variables are needed to be added to 
obtain a linear stochastic different equation that fit well to observations. There, 
the observed prices and the hidden variables can be considered as a vector bundle. 

There are many other practical examples where the measurements can not ex- 
actly be modeled. For instance in geophysical measurements implemented using 
explosives or air guns in water, the model for measurement process in not very ex- 
act. Thus the methodology of this paper, the combining of a set of measurements 
which are not well known, can be useful for improving reconstruction techniques 
with such inaccurately known measurements. 



AN INVERSE PROBLEM FOR A HYPERBOLIC SYSTEM ON A VECTOR BUNDLE 5 



The rest of the paper is organized as follows. In Section 2 we present some 
facts related to vector bundles, rigorously formulate our basic assumptions on 
sources used in energy measurements and state our main result. In Section 3 we 
explain that our assumptions on sources are quite natural by constructing such 
sources as sequences of realizations of suitable independent identically distributed 
Gaussian random variables. Section 4 is devoted to reconstruction of a bundle, 
the Riemannian structure on it and the system from measuring how much energy 
is needed to produce different combinations of the sources. 



2.1. Real vector bundles. Let M be a closed connected smooth m-dimcnsional 
manifold. Let V he a smooth real vector bundle over M and denote hyn-.V^M 
the projection onto the base manifold. Each fiber 7i~^{x) has the structure of a 
real vector space isomorphic to R*^ and a Riemannian structure that assigns a 
smooth inner product {■,-)x- 

Given two vector bundles tti : Vi — > M and 7r2 : V2 — > M, a map $ : Vi — > ^2 is a 
vector bundle homomorphism if tti = 7r2$ and $ restricts on each fiber to a linear 
map 7rf^(a;) 7r^^(a;). A bundle homomorphism $ : Vi ^ V2 with an inverse 
which is also a bundle homomorphism is called a vector bundle isomorphism, 
and then Vi and V2 are said to be isomorphic vector bundles. If $ preserves the 
Riemannian structure, then $ is called an isometry. 

Let {{Ua, 0a)} be an atlas of a vector bundle V, i.e. {[/«} is an open cover of M 
and (pa '■ 7r~^(f/a) ^ f/„ x M*^ are local trivializations. For any two vector bundle 
charts {Ua,<t>a) and {Uf3,(l)/3) such that f/^ H f/^g 7^ 0, the composition (pa o (pp'^ '■ 
(UanU/s) xW^ ^ (UanU/s) xW^ must be of the form 0„ o0^^(a;, t;) = {x,ta0{x)v) 
for some smooth map t^/? '■ Ua^Uf^ — > GL((i, M) which is called a transition map. 
For a given vector bundle atlas the family of transition maps always satisfy the 
following properties: 



A family of maps {tajs], tafs : Ua^Up ^ GL(d, R) which satisfies (2.1) for some 
cover of M is called a GL{d,W)-cocycle. 

The following theorem gives the minimal information required to reconstruct a 
vector bundle up to an isomorphism. 

Theorem 2.1. [21, Section 9.2.2] Given a manifold M , a cover {U a} of M and 
a GL{d,M.)-cocyde {tais} for the cover, there exists a vector bundle with an atlas 



2. Definitions and main results 



taf){x) ot/3a{x) 
taj{x) O t^f){x) O tpa{x) 



Id for all X G Ua and all a; 

Id for all a; G f/„ n Uf^; 

Id for all X e UaHUpnUj 



(2.1) 
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{{Ua,<f)a)} satisfying (pa^ (t>j3^{x^v) — {x,tai3{x)v) on nonempty overlaps Ua^Up. 
A vector bundle constructed in such a way is unique up to an isomorphism. 

2.2. Function spaces. Denote by C°^{M,V) the space of smooth sections of 
the bundle V. The inner product in L'^{M, V) is defined by 

(</.,V^)= / {(t>{x)M^)),dV{x) 
Jm 

where dV{x) is a smooth measure on M. For the vector bundle V over M, we 
consider also a related vector bundle V over M x R. When V has the transition 
functions ta/3 : Ua Ug ^ GL{d,M.), the bundle V has the transition functions 

tap ■■ (Ua X n (f//3 X R) GL{d,R), ta,f3ix,t) = taf3{x), {x,t) E M X R. Let 
L'^{M X R, V) and H'^{M x R, V) be the LP' and Sobolev spaces of sections of the 
vector bundle V , respectively. Wc often write L'^{M xR,V) = L'^{M x R), etc. 

Now let U be an open set on M and / be an open interval in R. Then we set 

H%U X /) = {F|t/x/ : F G H'{M x R)}, s G R, 

H'{U x /) = {F G H'{M X R) : supp (F) C THH}, s G R. 

For the spaces L^(C/ x I), H^{U x /) and C°°{M x /), sometimes we identify 
functions with their zero continuations in M x R. We recall the duality 

{H'{U X I))' ^ H-%U X I). 

2.3. Formulation of the problem and main result. Let 

A:C^{M,V) ^C'^{M,V), yl = ^ a„(a;)a^, 

|q|<2 

be an elliptic self-adjoint positive second-order partial differential operator on sec- 
tions of real vector bundle over a closed smooth manifold. In local trivializations, 
aa{x) are matrices. Consider the hyperbolic initial value problem 

{d^ + A{x,d:,))u{x,t) = F{x,t), xeM, t>-T, ^2 2) 

U{X, t)\t=-T = 0, dtu{x, t)\t=-r = 0, 

with some large r > and supp (F) C M x (— r, -|-oo). Let us denote by u^{x, t) 
the solution to (2.2). 

The energy of the wave is defined by 

Ea{F, t)^ \ f {{dtu^ix, t), dtu^{x, t)), + {Au^{x, t),u^{x, t)),)dV{x) 

where dV is a smooth measure on M. Since the operator A is self-adjoint, the 
energy is conserved, that is, if supp (F) C M x [to^^ol ^' i^" > then 
EA{F,t')^EA{F,t"). 
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Let Bm — {Up,p — 1,2,...} be a basis of topology of M and assume that 
Ui — M. Let B]^ = {Iq, 5 = 1, 2, . . . } be a basis of topology on M consisting of 
open intervals. By basic experiments we mean producing the following sets of 
sources: 

(BE) for any Up G Bm and G Br, we produce a set of sources ^Up,iq = 
{F^'^, /c = 1, 2, . . . } supported on Up x Ig which is dense in C{Up x Ig), 
\iv+ \ or in C°°(M x 7g), if p= 1. 

In Section 3 we show that the basic experiments can be almost surely gener- 
ated by taking some sequences of realizations of suitable independent identically 
distributed Gaussian random variables. 

Denote by {TTF){x,t) = F{x,t + T) the time shift for a source F. Let J G N and 

^ ^UpU)'hu}^ ^pU) ^ = (^^(i)'Co)) e Bm, J = 1, . . . , J. Assume that 

for any <t e Q and < 7} G Q, we can measure energy 

J 

EAiJ2^T,Fj,t) 

j=i 

provided that the time intervals [t"^^ — Tf^f^^..^ — Tj] are disjoint. These mea- 
surements are called admissible measurements. This means that we can combine 
separate basic experiments, delay the sources in time, and measure how much 
energy is used in the combined experiments. 

In the paper using admissible energy measurements for the sources generated by 
the basic experiments, we give a method to recover a bundle V over the known 
manifold M (up to an isomorphism), the Riemannian structure {•,-)x on it and 
an operator A. As our reconstruction procedure consists of infinite number of 
steps we formulate our main result as an uniqueness result. 

To simplify the notations, for countable bases Bm and B^, we renumber the 
sequences of sources generated by the basic experiments using one index, i.e. 
Fj = Fl\i^''^^^\ j G N. Letting /sT C N be a finite set, we denote 

Thus, we have the following uniqueness result. 

Theorem 2.2. Assume that a closed smooth manifold M together with its topo- 
logical and differential structures and a smooth measure on M are given. Let V 

and V be realj^mooth vector bundles with Riemannian structures over M and 
A{x, dx) and A{x, d^) be elliptic self-adjoint positive second-order partial differ- 
ential operators on sections of V and V , respectively. Let {Fj)'jLi = (-f^'m 



8 



KRUPCHYK AND LASSAS 



and = {F^lj^'''^-'^)'^-^^ be sequences of sources generated by the basic exper- 

iments for V and V , correspondingly. If for any i^T C N, (T,)^^, < Tj G Q 
andQ<teQ, 

Ea(K, (Tj)^=„t) = Ej(K, (7;)~ „ t), 
then there is an isometry ^ : V ^ V such that A — 

Roughly speaking, Theorem 2.2 states that if we do measurements using combi- 
nations (i.e., sums) of some unknown basic sources and find the energies 
needed to produce all combinations of these sources, then the physical model 
(i.e., the Riemannian bundle structure and the operator) is uniquely determined 
up to an isometry. 

3. Generation of basic experiments. Random sources 

The purpose of this section is to prove that for any 7 G Br and U G Bm, 
sequences of sources, which are dense in C°°{M x /) and C{U x I), can be almost 
surely generated by taking some sequences of realizations of suitable independent 
identically distributed Gaussian random variables. This shows that the sequences 
of sources we consider are generic. 

Let {fl, E, P) be a complete probability space. Recall the notion of a Gaussian 
random variable taking values in a real separable Hilbert space Y with inner 
product (• , ■ )y- Let us identify the dual of Y with Y using Riesz representation 
theorem. Let A/y be the Borel cx-algebra of (F, r"") with r"" being the weak 
topology of Y. Note that the separability of Y verifies that A/y coincides with 
the Borel cr-algebra of (F, t""), where r" is the norm topology of Y. An Y -valued 
random variable F is a measurable map F : (Jl, E) — > (Y, A/y)- In this paper we 
consider only random variables with values in separable Hilbert spaces. We say 
that F is Gaussian, if for any ip &Y the real- valued functions u ^ {ip, F{uj))y are 
Gaussian random variables. Then F-valued random variable has the expectation 
EF G y and the covariance operator Cp — Cf,y :Y ^Y {m. the space Y) which 
are defined by 

(EF,V')y = E(F,V')y, for ^ e 

and 

E ((F - E F, 0)y (F - E F, V')y) = (Cj^0, V>y , for (t>,^ eY. (3.1) 

Note that if Yi and Y2 are two separable Hilbert spaces so that Yi C ^2, any Yi- 
valued gaussian random variable F can be also considered as l2-valued random 
variable but the convariance operator of F in the space Yi, Cp^i, do not generally 
coincide with the convariance operator of F in the space Y2, Cf^y2j but 

Cf,Y2 — JCf,YiJ* 



AN INVERSE PROBLEM FOR A HYPERBOLIC SYSTEM ON A VECTOR BUNDLE 9 

where J : Yi ^ Y2 is the identical embedding and J* : Y2 ^ Yi is its adjoint. 
The covariance operator Cp : Y ^ Y is always a selfadjoint non-negative trace 
class operator inY [4, Theorem 2.3.1]. 

Let F be a Gaussian random variable in Y and let Cp '■ Y ^ Y have the 
eigenvalues Xj > and let ipj, j = 1,2, .. . be the orthonormal eigenf unctions 

which form a basis in Y. Then F can be writen in the form 

00 

i=i 

where Gj = {F,(pj)Y, Gj : fl —>■ M. are independent normalized Gaussian random 
variables and aj = ^/Xj- 

Lemma 3.1. Let Fj, j = 1,2, ... , be independent identically distributed Gaussian 
Y -valued random variables with EFj = and covariance operators Gp. : Y ^ Y 
being injective. Then the set {Fj{u!),j — 1,2, . . .} is almost surely dense in Y 
with weak topology. 

Proof. To see the claim one has to show that for any open set 7^ in weak 
topology, 

P(u~i(F,- e W)) = 1. 

Let Tz^r,e ^ {y eY : \{z, y)Y - r| < s} for 2; e y, r e IR and e e R+. Then the 
set 

is a basis of [Y, r^) . 

Let us consider an arbitrary set A e Bbasis \ {0}, i-©- 

N 

^ = r\'^z„n,e„ ZiEY, neR, SieR+, NeN (3.2) 
1=1 

Let us next consider a single random variable Fj. Our first goal is to prove that 

F{Fj e A) > 0. (3.3) 

Let P : y ^ y be an orthogonal projector onto a finite dimensional space 
Z = span (zi, . . . , Zi\f) and endow Z with the same inner product as Y. 

Since the covariance operators Cj?. : F — > F is injective, Gp^ : Y Y has the 
following positivity property: 

iftpeY satisfies {Gp^i/j, = 0, then ^ = 0. (3.4) 

Then PFj is Z-valued random variable having zero expectation and covariance 
operator 

CpFj ■ Z ^ Z, Gpp. = PGp.P. 
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Now a t/j E Z is such that {Cpp.t/j, t/j)z — 0, then 

As Pip = ip, we have {Cfj1P,iP)y = and thus ip = 0. This imphes also CpF^ 
Z ^ Z satisfies the positivity property (3.4). 

Now as A is open and non-empty, PA C Z is open and non-empty. Identifying 
Z with R^, (3.4) yields that Cpp^ is an invertible matrix and thus, we have 

FiPFj e PA) I exp(-^x • Cpl,x)dx > 0. 

JpA 2 ^ 

By definition of P and A we get that Fj e A if and only if PFj G PA. Hence, 

P(Fj eA)= F{PFj e PA) > 
that shows the claim (3.3). 

Now let be a dense set in (F, r") with norm topology and consider the 

enumerable set 

5~ ^ {n^^^T,,,,,,^ : y, e {y^}r=l,r^ eQ,e,E Q+,Ne N}. 

Denote by Bn = {y e Y : ||y||y < R} the ball of the radius RinY and introduce 
enumerable sets 

Vi ={V nBp: ReQ+,V e 5™}, 
V3 ^{v e V2 : r{Fj e V) > 0}. 

Now let A be a set with the following property 

AeVi and A G A. (3.5) 

Let us number all the sets with the property (3.5) by Ai, A2, Then 

00 

We set 

n 

Vn^[jAi, 

i=l 

and thus, Ki C Ki+i, Ki £ V2 and Vn C ^1. By the monotone convergence 
theorem, 

lim P (F. e 14) = P (^j e A) > 0. 

n— >oo 

Hence, for a non-empty set A of type (3.2), there exists a set F e V3 such that 
V dA. 
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By Kolmogorov's 0-1 law, we have for any V E 

P(u;^i(F, G V)) = 1, i.e. F{[\Jf^^{Fj e V)Y) = 0. 
Since the set V3 is enumerable, we get 

VeVa 

Thus, we have 

1 = nnvevs u~ 1 {Fj e V)) = P({Vy eVs3jeN: Fj e V}). 

As any non-empty weakly open set 1^ in y contains a non-empty set A of the 
type (3.2), and any such set A contains a set F e V3, this yields 

F({yW er'",Wy^ 0, 3j e N : Fj e W}) = 1. 

Thus the set {Fj : j = 1, 2, . . . } is almost surely weakly dense in Y. 

□ 

As compact operators map weakly convergent sequences to norm convergent se- 
quences, we obtain the following corollary. 

Corollary 3.2. Let Fj, j G N, be independent identically distributed Gaussian 
Y -valued random, variables with KFj — and covariance operators Cp^ : Y ^ Y 
being injective. Let Y be a Hilbert space so that Y Y is compact and Y G Y 
is dense. Then the set {Fj{u!),j = 1, 2, . . . } is almost surely dense in Y in norm 
topology. 

Next we give an example how one can construct a set of the Gaussian random 
variables which is weakly dense in H'^{M x M) for any s > 0. First recall that 
the Gaussian random variable is determined uniquely by its expectation and 
covariance operator. An operator Cf,h'> ■ H^{M x M) — > H^{M x M) is a co- 
variance operator of some random variable F if and only if Cf,h^ is self-adjoint 
non- negative trace class operator in H^{M x M), see [4, Theorem 2.3.1]. 

Example 3.3. Let us construct the Gaussian random variable F which can be 

considered as a random variable in the Sobolev space H^{M xM.) with an arbitrary 
s > such that KF = and it has an injective covariance operator. First assume 
that s = and consider a Schodinger operator 

L = {-Ag - dl + 1^) : L^{M x M) ^ L^{M x M). 

It is known [9] that the spectrum of L is discrete. We set Cf,ho — ■ 

Now assume that s > and let G be a Gaussian random variable in H^{M x M) 
with the covariance operator Cg,H' : i?*(M x E) ^ H^{M x E) and EG = 0. 
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Let J : H^{M X R) ^ L^{M X R), J(f) ^ (f), be the embedding. Then the random 
variable JG : (Q, E) — > L^{M x R) has the covariance operator 

CjG,L^ = JCcHxJ*- (3.6) 

One can see that J* = K'' : L'^{MxR) ^ H'{MxR) where K = {-A-d'^ + 1). 
Now setting Cjg,l^ = e'^ , (3.6) implies that Cg,h^ = e'^K'^ : H'{M x R) ^ 
H^{M X M) which is clearly injective. 



Let Fj, j = 1,2, . . . be independent identically distributed Gaussian H'^{M x M)- 
valued random variables, s > 0, with EFj = and injective covariance operators 
Cf^ : H'{M xR)^ H'{M x R). Then Lemma 3.1 the set {Fj{u),j = 1,2,...} 
is almost surely weakly dense in H^{M x R) for all s > 0. Let / e Br. Then 
{Fj{uj)\Mxi: j — 1)2,...} is almost surely weakly dense in H^[M x I) for all 
s > 0. Consider sq, s > sq > 0. Then by Corollary 3.2, {Fj{Lj)\Mxi, j = 1, 2, ... } 
is almost surely dense in H^^^M x I) with norm topology for any sq > 0, and 
thus in C°°{M x I) with a Prechet topology. Recall that the space C°°{M x 1) 
is a Frechet space with the topology determined by the family of seminorms 
Pk{F) = ||F||cfc(Mx7)> A; = 0, 1, 2, ... . 

Let U e Bm- Then 

{Fj{u;)\u.iJ = h2,...} (3.7) 

is almost surely dense in C(C/ x 7). Also, the set (3.7) is almost surely dense in 
L^(C/ X 7) with norm topology. 



4. Proof of Theorem 2.2. 

4.1. Measurements. For convenience of the presentation, we use the following 
notations. Let Ii = {ti,tt) and I2 = (^2 '^2) intervals in M and Tq e M. We 
write 

7i < I2 if tf < t2, 

7i < To if tf < Tq. 

In this section we discuss what kind of information can be computed using ad- 
missible measurements. For this end, we need the following lemma. 

Lemma 4.1. (i) Let Fi G C{Ui x U), Ui G Bm, L, G B^, i = 1,2 (recall that 
we identify Fj with their zero continuations) . Then for any time Tq, 

Ea{Fi ± F2, To) X llii^^To) ± u^'{to)\\%i^m,v) + Ptu'^' (to) ± dtU^'{ro)\\h^M,vy 

(4.1) 
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(ii) Let Fi e C°°{M x 7i), U e Bm, i ^ 1, 2. Then for any To, h < tq, i ^ 1, 2, 
and s — 1,2, ... , 

(4.2) 

where fi x /2 if there are constants ci, C2 > such that ci/i < /2 < C2/1. 

Proof. The definition of the energy immediately impHes (4.1). Let us prove (4.2). 
As Ii < To, we have dfU^^{To) — —Au^'^{to) and moreover, 

{dt)V{To) = {-l)'AV^{To), k=l,2,.... (4.3) 



As A is an eUiptic second order operator and A = is not its eigenvalue, we get 

ll^'^'^llL2(M,y) ^ Il'^llff2fc(jv^^y), k — 1,2, ... , 



k k k 2 2 (^•^) 
{A{A u),Au)>^ \\A U\\jjl(^f^jy^ X ||'w||^2fc+l(jy^y), k — 0,1,2, 

Now for s = 2A; + 1, /e = 0, 1, . . . , (4.3) and (4.4) imply 
EA{{dty{Fi-F2),To) 

= ^{{A{A\dtU^^iTo) - dtU^^iTo))),A'{dtU^^{To) - dtU^'iTo))) 

+ \\A'-'\u^^{To)-u^'^{romhiM,v)) 

For 2k, k^ 0,1,..., (4.3) and (4.4) yield 
EA{{dty{F,-F2),To) 

= 1{{A{A\u^^{to)-u^^{to))),A\u^^{to)-u^%to))) 

+ \\A%{u^^{To)-U^%TomhiM,V)) 

X \\u^'{to) - 1i^'(To)||1^2fc+i(M,y) + \\dtU^'{To) - dtU^'{To)fH2k^M,vy 

□ 

In a sequel to prove that convergence of sequences of sources Fj — > F when j — > oo 
implies convergence of waves u^^{To) u^{Tq) and dtU^^lTo) dtU^iTo) when 
J — > oo in appropriate spaces we need the following estimates. 

Lemma 4.2. (i) Let either F e C~(Mx7o) andO < s or F e C(Ux7q) 
and s = 0. Then for any Tq, 



U iTo)\\Hs+i(M,V) < C'/o,To 11-^11 L2(7,j^^s+i 



\\dtU^{To)\\H^(M,V) < C/o,ro||i^||L2(7o,//^ 



(4.5) 



{M,V))- 
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(ii) Let F e L'^(Iq, H\M, V)), Iq = {to, t^). Then the function t Ea{F, t) 
is continuous on any finite time interval Ii — (tQjtf). 

Proof. The proof of (i) is quite well-known for hyperbolic systems but we include 
it for convenience of the reader. As A is an eUiptic self-adjoint operator, there is 
an orthonormal basis of eigensections {</?;} of A in L^(M, y). Denote by A; > an 
eigenvalue of A which corresponds to (pi. Thus, the wave can be represented 
by the Fourier series 

oo 

u^{x,t) = J2^^it)M^)^ <W = {^^{■,t),ipi) (4.6) 
1=1 

with coefficients satisfying 

ufit) = £ F,is)ds, F,is) = (F(, .), ^,). (4.7) 



Thus, 

\uf{To)\'<Cj,,T,f^_ \Fi{s)\'ds. 



Hence, we have 

oo 



1=1 

„f+ oo 



•^*o 1=1 



— Cio,To l|-^("S)||H«+i(M,y)C^"S — C'/o,Toll-^llz,2(7g_j^s+i(M,y))- 

The second estimate in (4.5) can be proven in the same way. Now uniform 
convergence of series (4.6) and (4.7) implies (ii). □ 

RemEirk 4.3. Lemma 4-2 implies that measuring energy EA{^j^x'^TjFj,t) for 
Tj e Q and t & Q, we can find the energy EAC^j^j^TT^Fj, t) for any Tj e R and 
teR. 



The lemma below shows that for any source F defined for a time interval Iq, there 
exist sequences of sources generated in the basic experiments and defined for a 
time interval Ii satisfying Iq < Ii such that the sequences of the corresponding 
waves approximate the waves and at time Tq > /i. 

Lemma 4.4. Let either F e C°°{M x /q) and < s e R or F e C{U x Jq) 
and s = 0. Then for F and any tq, Iq < tq, there are sequences of sources 
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Fj, Gj e TM,h, h e -Br, j = 1, 2, . . . , 7o < -^1 < To, such that 

\im{u^^{To),dtU^^{To)) = {u-^{To),dtU-^{T^)) lu H'+\M,V) X H'{M,V). 

j— >oo 

(4.8) 

lim («^^(to), dtu'^^iro)) = {u^{to), dtU^{To)) m H'+\M, V) x H'{M, V). 

(4.9) 

Proof. Consider first the case when F G C'^{M x /q). Let e C°°(M) be such 
a function that t/j{t) — 0, when t < ti, and 'ip{t) — 1, when t > t^. Then 
= —t/ju^ is a solution to (2.2) with the right-hand side F — —^dtu^dti^ — 
u^dli) e C^{M X R),_supp (F) c M X 7i. Thus, F e C°°(M x Ji) if the 
functions from C°°(Mx Ji) are identified with their zero continuations. Moreover, 
{u^ ijo) , dtu^ {tq)) = {u^^{To),dtU~^{To)). By (BE), there are sources Fj G J-'mji 
such that Fj F in C°°{M x 7i) when j oo. Hence, (4.5) implies (4.8). The 
existence of a sequence Gj G Tmji which satisfies (4.9) can be seen in the same 
way. The case when F E C{U x Iq) is proven similarly. □ 

Lemma 4.5. Using admissible measurements, for a given source F G J^ujo or 
F G J-M,io, one can determine sequences of sources Fj, Gj G J-M,h, j — 1, 2, ... , 
which satisfy Lemma 4-4- 

Proof. By Lemma 4.4 there exists a sequence of sources Fj G J-'mji, J = 1, 2, . . . , 
which satisfies (4.8). Thus, (4.1) implies that 

lim EA{Fj + F, To) = 0. (4.10) 

j— ►oo 

Now since by our assumption energy (4.10) can be evaluated using admissible 
measurements, we can verify if an arbitrary sequence of sources Fj G J-'mji 
satisfies (4.10) or not. Thus, we can choose a sequence of sources Fj G TmJi 
which satisfies (4.10). Hence, (4.1) yields that (4.8) holds for this sequence. 

By Lemma 4.6 below, the sequence Gj G J-M,h which satisfies (4.9) can be 
determined in the same way. 

□ 

Recall that by our assumption we can only compute energy of the sums of sources 
generated in the basic experiments. The following lemma shows that we can 
compute the energy of the difference of the sources. 

Lemma 4.6. Using admissible measurements, for any sources Fi G Tij^^i-, i = 
1, 2, the energy Ea{Fi — F2, tq) can be computed at any time Tq, li < Tq, i — 1, 2, 
and also at any time Tq G I2 provided Ii < l2- 
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Proof. Let us prove the first part. Let I2 G be such that li < I2 < tq, i = 1,2. 
Then by Lemma 4.5, one can determine a sequence of sources F? e ^mT2 ^^^^ 

hm («^/(to), dtu^'iro)) = («-^^(ro), dtU-^'{To)) in H\M, V) x L2(M, V). 

Thus, 

- F2, To) - hm £;^(Fi + Ff, To). 

Since the supports of and Fj are disjoint in time, the above energy can be 
computed using admissible measurements. 

Assume now that Ji < I2 and tq G l2- Then as in the first part of the proof, 
fixing some arbitrary interval Ii G B^, such that Ii < Ii < I2, one can determine 
a sequence of sources Fj G J^j^ safisfying 

lim E^(F/ + Fi,To) =0. 

Thus, one can compute 

Ea{F, - F2, To) = lim Ea{-FI - F2, To) = lim Ea{F} + F2, Tq). 

J— »oo J— >oo 

□ 

Lemma 4.7. For any sources Fi G ^Uiji, i — 1,2, T > and tq, Ii < tq, 

i = 1,2, the energy 

Ea{Fi - TtF2,Tq) 

can be evaluated using admissible measurements. 

Proof. Let I3 G be such that Ii < Is < Tq, i = 1,2. Then a sequence of sources 
Ff G J^M,i3 satisfying 

lim F^(F? + ttF2,To) = 

j-*oo 

can be determined in the same way as in Lemma 4.5. Now using admissible 
measurements, we can compute 

Ea{Fi - ttF2, To) = lim Ea{Fi + F/, To). 

□ 

Corollary 4.8. For any sources Fi G J^Ui,ii, Ti > 0, Qi e Q, and Tq, Ii < Tq, 
i — l,...,p, the energy EA{Yli^iqiTTiFi,To) can be computed using admissible 
measurements. 
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Proof. Assume for simplicity that Qi — rrii/ni > and I — n^^^rij. Then 
EA{Y7i=i(li'^TiFi,TQ) = l-'^EA{J2i^rhiTTiFi,To) with fhi = rriiUk^ink G N. Now 
setting TrjpkF^ — TTiF^, k = 1, . . . , mj, we write 

P P fhi 

Ea{Y1 fhiTTiFi, To) = Ea{J2 Y1 ^T^Fi, To). 
i=l i=l k=l 

Taking arbitrary intervals if G B^. such that 

/i </!<••• < /r <ii< ...if' <■■■ <i^' < To, 

and using admissible measurements, one can find sequences of sources Fj *^ e 
Tj^ ik such that 

lim£;^(Ff -T^.^;^To) = 0. 

Thus, 

p V fhi 

EAiY.'ii^TA^ro) = hm F^(^^Ff ,ro). 



j=l i=l ik=l 



□ 



Corollary 4.9. For any sources Fi G J^u^ja Ti > 0, hi & M, and tq, Li < tq, 
i = l,...,p, the energy EA{Y^^=ihiTTiFi,To) can be computed using admissible 
measurements. 

Proof. Taking sequence /i" G Q such that /i" hi when n — > oo, we can compute 

p p 
EA(ThiTTiFi,To) = lim V /i^T.i^i, To). 

□ 



n— ►oo 

1=1 1=1 



Lemma 4.10. For any source F G J^ujo CLnd time tq, Iq < tq, the energy 
EA{d^F, Tq) of its s-th time derivative at time tq, s — 1,2, ... , can be computed 
using admissible energy measurements. 

Proof. Using the forward difference approximation to the s-th order derivative 
and Corollary 4.8, we can compute 

EAidtF,To) = Ihn -^£;^(^(-l)'= (l) ^s-k)hF,To) 

~^ k=0 ^ ^ 

where are binomial coefficients which are natural numbers for all s and k. 

□ 
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Lemma 4.11. For any sources F,H& J-u,io o-i^d time Tq, Iq < Tq, the energy 
EA{d^{F — H),To), s — 0,1,2,..., can be evaluated using admissible measure- 
ments. 

Proof. The proof follows from the forward difference approximation to the s-th 
order derivative and Corollary 4.8. 

□ 

Using the forward difference approximation to the s-th order derivative and Corol- 
lary 4.9, we have the following lemma. 

Lemma 4.12. Let Fi e ^Ui,ii, i = 1, 2. Then for any s = 1, 2, . . . , the energy 

EA{d^,F,-F2,To) 

at any time tq, li < tq, can be computed using admissible measurements. 

The following lemma shows that admissible measurements allow us to compute 
inner products between sources generated in the basic experiments and the time 
derivatives of waves corresponding to such sources. 

Lemma 4.13. Let Fi e J-M,h '^"'^ -^2 G < h < h- Then for any 

I'd G h, the inner product 

{F2{ro),dtU^^iro)) 
can be computed using admissible measurements. 

Proof. Since A is self-adjoint operator, we have 

dtEA{F, t)= [ {{dlu^{x, t), dtu'^ix, t)), + {Au^{x, t), dtu'^ix, t)),)dV{x) 

J M 

= / {Au^{x,t) + dlu^{x,t),dtU^{x,t))^dV{x) 

= [ {F{x,t),dtU^{x,t)),dV{x). 
Jm 

Thus, we get 

dt[EA{F, + F2, t) - Ea{Fi - F2, t)] = 

2 / {{F,{t),dtU^^{t)), + {F2{t),dtu'''{t)),)dV{x). 
Jm 

As Ii < I2, we have 

supp {u^''') n supp (Fi) = 0. 

Hence, 

dt[EAiFi + F2,t) - Ea{Fi - F2,t)] = 2{F2it),dtu'^'{t)). 
Thus, by our assumption and Lemma 4.6 the inner product {F2{to) , dtu^^ {tq)) 
can be evaluated using admissible measurements. 
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□ 

Corollary 4.14. Let Fi e J^M,h and F2 G ^u,i2> —T<Ii<l2- Then for any 
To e I2 and Tq, I2 <To, the inner product 

{F2iTo),dtu''^{To)) 

can be computed using admissible measurements. 

Proof. The claim foUows by applying the proof of Lemma 4.13 with ttq-toFi 
instead of Fi and the fact that M'^^o-^o^i(ro) = tto-toU^^{tq) = u^^{Tq). The 
latter is a consequence of (4.7). 

□ 

The lemma below allows us to compute inner products between sources generated 
in the basic experiments and waves corresponding to such sources using admissible 
measurements. 

Lemma 4.15. Let Fi G J-'m.Ii and F2 G J-u,i2) < h < h- Then for any 
Tq G I2 and To, I2 < Tq, the inner product 

{F2{to),U^^{To)) 

can be computed using admissible measurements. 

Proof. Let Is G B^. be such that Ly < I^, < h- Then for Fi G Tmji, there is a 
sequence of sources Gj G J-m,Iz satisfying 

lim£;^(9tG',-Fi,To)=0. (4.11) 

Since Ji < Is < To, by Lemma 4.12 the energies in (4.11) can be computed using 
admissible measurements and hence, we can verify for any sequence (G'j)j?,i, 
Gj G Tmjz, if (4.11) is valid or not. Let take the sequence Gj which satisfies 
(4.11). Thus, (4.11) and (4.1) imply that 

lim (atK^^(To),a2w«^(To)) = {u^'{To),dtu'^'{To)) in H\M,E) x L\M,E). 
By Corollary 4.14 we can compute the inner products 

(F2(To),at«^^(To)). 

Thus, 

{F2{TQ)y\%)) = lim(F2(To),9tti«^(To)). 

3^00 

□ 

To reconstruct an operator A we need the following corollary which allows us to 
evaluate inner products between sources generated in the basic experiments and 
any order of time derivative of waves corresponding such sources. 
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Corollary 4.16. Let Fi e J^mji o-n-d -F2 £ ^u,i2! —T<h<l2- Then for any 
To e I2, Tq, I2 < Tq and s e N the inner product 

can be computed using admissible measurements. 

4.2. Generalized sources. Let Iq e Br- For any s = 0, 1, 2, . . . , let us introduce 
the following space of sequences of sources 

n,rc = {(^.)r=i ■ e ^M,/o, .lim EAidtiF, - F,),To) = 0}, Jo < Tq. 

j,fe— >oo 

By Lemma 4.11, the space ^f^^ro ^® constructed using admissible measure- 
ments. 

The definition of the space J^j^^ j-o (^-2) imply that {u^^ (tq) , dtu^^ {tq)^°°__^ is a 
Cauchy sequence in H'+\M,V) x H'{M,V). Since H'+\M,V) x H'{M,V) is 
complete, there are a G H^^^{M, V) and h E H^{M, V) such that 

lim u^^{to) = a in H'+\M, V), 

j-»oo 

lim dtu^^ (to) = 6 in i/*(M, V) . 

Thus, we can define an operator 

W : ^ H'+\M, V) X i/^(M, \^), 

= (lim M^^(ro), lim 9,;/^(ro)). 

J— >oo J— >oo 

Moreover, we can define a semi-norm on the space ^f^^Tg of sequences of sources 
Lemma 4.10 allows us to compute the semi-norm (4.12). 

We say that (-Fj)^i G J^i^^ro 1^ equivalent to (-Fj)jLi G .^/q^t-q denote by 
(F,)- , ~ (F,)- , if 

lim u^^ (to) = lim u^^ (to) in (M, T/) , 

J— >oo 

lim atM^^ (To) = lim dtU^^{To) in i/^(M, 1/). 

J— >oo J— >oo 

This is equivalent to that fact that \\{Fj — Fj)'^^-y\\jrs^ =0. Further, we define the 

space J^i^^rol ~ (4-12) becomes a norm on this space. Finally, we complete 

^iQ^ro/ ~ ^1^1^ respect to the norm (4.12), d{!Ff^^^/ ~) = ^Io,to- "^^^^ space is 
called the space of generalized sources. 

Remcirk 4.17. The space ^^j^^ro of generalized sources can be constructed using 
admissible energy measurements. 
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By the definition of the completion, the space ^/q^t-q consists of the sequences 
F — {{Fj)^^)'^^-^ which are the Cauchy sequences with respect to the norm 
(4.12). We define 

fc— >oo J— >oo 

dtU^{To) = hm hm dtU^^{TQ), 

W{F) = {J{To),dtJ{To)). 

By the definition of the space T^j^^^^ the map W : ..^ H'+^{M, V) x 
i7^(M, is injective. 

Lemma 4.18. W : T]^.^^ H'+^{M,V) x H'{M,V) is surjective. 

Proof. First we show that for any (a, b) G C°°(Af, V) x C'^{M, V), there is F G 
C°°(M X Jo) such that u^{tq) = a and dtu^ {tq) = b where is a solution to 
(2.2). Indeed, let w be a solution to 

(dt +A)v = 0, in M X M 

v\t=To = Ct, 
dtV\t=ro = b. 

Let tp G C°°(M) and ipit) = ii t < and ipit) = 1 if t > t^. Then u^{x,t) = 
v{x,t)ij{t) is a solution to (2.2) with F = 2dtvdti) + vd'^^) G C°°{M x R), 
supp (F) C M X To. Thus, F G C~(M x Jq), if the functions of C~(M x Jq) are 
identified with their zero continuations. Moreover, u^{to) = a and dtu^ijo) = b. 

Now let (a, b) G H'+\M, V) x i:f^(M, V). Since C°°(M, V) is dense in all Sobolev 
spaces, there are Oj, bj G C°°{M, V) such that 

a = lim in H'+^{M,V), 

6= lim bj in //"(M,!/). 

By the first part of the proof, there are sources Fj G C°°(M x /q), such that 
u^^i^o) = cij and dtU^^{To) = bj. Since ^m,/o is dense in C°°{M x /q), for any 
there is F^ G .Fm,7o such that limfc^ooF^^' = F^ in C°°(M x Jq). Thus, (4.5) 
yields that limfe_^ooii^fc(To) = % in H^'^^(M,V) and limfc_,oo 9tii-'^fc(ro) = bj in 
-/^^(M, "K). Hence, lim^^fc^^oo '^^^^(to) = a in H^~^^(M,V) and limj^fc^^oo •^f'U'^fe (tq) = 
b in H'{M,V). Moreover, (4.2) implies that limfc,i_,oo ^^^(^t (F^ - F/),ro) = 
and hence, (F^)- , G J^^,^. 

□ 
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Now we define the space 



T 

■J \ 



■oo 



n 



-fo>TO' 



■s 



s=0,l,2,... 




4.3. Bundle reconstruction. Let denote by S^^ the delta-distribution at Xq e 
M, i.e., 



The main point of the bundle reconstruction is to use admissible energy mea- 
surements to find for any point xq G M, sequences of sources that converge to 
delta-distributions XS^^, A e 7r"^(a;o). Recall that we denote by m the dimension 
of the manifold M and then 



S,,eH-"'/^-%M,V), ^ H-^/^-'{M,V), ee(0,l], |a| > 1. 



If m is even, we set Sq = m/2 G N and otherwise, Sq = m/2 — 1/2 G N U {0}. 

Since M is a compact differentiable manifold, it admits a Ricmannian metric. 
Let us denote by B{xo,r) an open ball on M at a point Xq G M with radius r 
with respect to this metric. Then B{xo, r) is an open set with respect to topology 
on M. Let (rj)^^, > 0, be a sequence such that ^ when i ^ oo. Then 
for any point xq G M, there is a sequence of {U^°)^i of the sets G Bm such 
that 



For any i E N and any interval Jq G B^, there are sources J^i,xo,io = {^i,k '■ k = 
1,2,...} which are dense in C{U," x Jq). 



with some section A. Here Fj^fc(j) G Ti^xojo- 

Lemma 4.19. can he constructed using admissible energy measurements. 

Proof. Let an interval Ii G B^ and Tq G M be such that — r < /i < /q < Tq. Then 
by Lemma 4.15 we can compute the inner product 




Xq G f/f° C B{xo,ri). 



Now let be a set of all sequences such that k{i) G M, G /oflQ 

and 



limF,,fc(,)(t,) = A4o in H-''-\M,V) 



{Fi,k{ti),u'"{To)) 
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for any Fi^k e ^i,XQ,io-, Ue IoDQ and F e ^Ti,to- Let us show that ti)iGN e 
if and only if there exists 

lim{Fi,k(i){ti),u^{To)) (4.13) 



for any F e ^°i[,to- Indeed, if the hmit (4.13) exists, then as W : ^/"^^-o 
if*°+^(M, y) X H^°{M,V) is surjective, the Banach-Steinhaus theorem imphes 
that 

a = hm Fi k(i){ti) exists in H-">-^{M, V). 

i— >oo ' 

Since supp {Fi^k(i)i'ti)) C ?7f° C B{xQ,ri), ^ when i —>■ oo and a is a distri- 
bution, we have that supp (a) = {xq} and thus, a is a finite hnear combination 
of the delta-distributions at Xq and its derivatives with coefficients from 7i~^{xo), 
see [12, ex. 5.1.2]. By the choice of Sq, the space H~^°^^{M,V) contains delta- 
distributions while not their derivatives. Thus, we have a — XS^q- 

Hence, to construct the set ^^^o determine all sequences {k{i),tij.^^ such 

that the limits (4.13) exist. □ 

Lemma 4.20. For any Xq e M, C^^^j^ ^ 0- 

Proof. Let x — > A(x) be a section on Ui°. Then consider the following sequence 
of functions 

Giix,t) = Xaau:^{x)xioit) e L'iur x 7o), 

where 



^ vol(C/r)= / IdVix), i^l,2, 



As the set of sources J-'i^xo,io produced in the basic experiments is dense in 
L'^{U^° X /o), for any e > and i G N there is k{i) G N, k{i) oo when 
i — > oo, such that 



' Jlo 

Thus, by Fubini's theorem there exists ti G Iq such that Fi^k(i){ti) ^^'^ 

\\Fi,k(i){ti) - Gi{ti)\\L2^u:°) < e. (4.14) 

As Fi^kit) e C(To X C/^°), one can find E lodQso that (4.14) holds. Hence, 

lim Fi^k{i){ti) = lim AaiXt/-o (x) = A^^^ in H-'°-\M, V) 
as i7^o+i(M, V) ^ C{M, V) by Sobolev's embedding theorem. 

□ 
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(2) 

Let be a sufficiently small neighborhood of the point Xq in M. Let be 
a set of all sequences ((k(x,i),ti)iQ^) such that 

1. for any x G f/rg, (A;(a;, i), tj)^^^^ G -^i^/g, i-e. G /q H Q and there are 

FiXx,i) e J^i,x,/o such that limi^oo Fi^k{x,i){ti) = A(a;)(5a; in H~^°~^{M,V). 

2. X — s> A(a;) is a C°°-smooth section in f/r;,. 

Lemma 4.21. For any xq G M and sufficiently small neighborhood Uxq of xq, 
the set Cjj^ ^ and C\j^ can be constructed using admissible energy mea- 
surements. 

Proof. The fact that lo ^ for ^ sufficiently small neighborhood of Xo 
is clear. 

Now it is known that A is a C°°-smooth section in Uxq if and only if functions 

: ^ R, K{x) = {X{x), (t>{x))x 

are C°° smooth for all G C°°(M, F). Let an interval h e -Sr and Tq G R be 
such that /i < /o < tq. Since : F^,.^ ^ C°°(M, x C°°(M, F) is bijective, 
to check whether the functions K(f, are C°° smooth for all G C°°{M,V) is 
equivalent to check that the functions 

{\{x),u^{x,Tq))x = {\{x)5x,u^{tq)) = lim(Fife(^j)(ij),ii^(To)) 

I— >oo 

are C°°-smooth for any generalized source F G . The latter can be checked 
using admissible measurements. 

□ 

Let us denote by ^ul^ the set consisting of all sequences {{ki{x, i),tl)i^f^) ^^^^ , 
{(kd(x, i),tf)ien)x^u^^ such that 

1. {{ki{x,i),t\)i^fi)^^j^^^ G C^^^j^, I = l,2,...,d, i.e. there are Fi^ki(x,i) e 

^i,x,io such that limj^oo = X''{x)dx with C°° smooth section A', 

l^\,2,...,d; 

2. A^(xo), . . . , A''(xo) G TT ^(^o) are linear independent; 

3. d is a maximum number such that the above properties can be met. 

Lemma 4.22. For any G M and sufficiently small neighborhood Uxq of xq, 
the set >C[^j^ Iq ^ ^'^'^ ^uIq Io ^^'^ constructed using energy measurements. 

Proof. The fact that jC^^ 7^ for a sufficiently small neighborhood Uxq of xo 
is clear. 
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Now it is obvious that X^{xo), . . . , A'^(xo) G 7r~^{xo) are linearly independent if 
and only if the functionals 

are linearly independent. 

Let an interval Ii G B^. and Tq be such that — r < /i < /q < Tq. Then by the 
construction of the generalized sources 

7r~^(xo) = {p{xo) : p : M ^ V is a. C°° smooth section} 

= KKro):Fe7=~J. 

Now the linear independence of Kxi, I — 1, ... is equivalent to linear indepen- 
dence of the functionals 



Ki : ^ R, Ki{F) = {\\x^)y{x^, To)),, = {X\xo)d,„ u^{to)) 

= lim(Fi,fc,(^o,i)(t-),'u^(ro)), l = l,2,...,d. 

I— >oo 

Thus, using admissible measurements one can check whether X^{xo), . . . , X'''{xo) 
are linearly independent. 

□ 

Lemma 4.23. Admissible energy measurements make it possible to reconstruct 

local trivializations : T^~^{Ua) — >■ [/^ x M*^ of the vector bundle V and the 
GL{d,M.)-cocycle {ta/s}, tap : Ua^^Up ^ GL{d,M.) between any two local trivial- 
izations. 

Proof. Let Xq E M he an arbitrary point and choose a neighborhood C M of 
Xq such that JC^^^^ j^ ^ 0, i.e. 

\\x)^, = \v^Fi^ki{x,iM). {i.H^^i)A)ien)x^u,^ e ^ul,M^ l = l,2,...,d. 

As sections A^, . . . , A*^ are C°° smooth, there is a neighborhood U C Ux, of xq 
such that {A^, . . . , A'^} is a frame field over U. 

Let Ii e Br and Tq be such that —t<Ii<Io<To. Now using admissible energy 
measurements, one can evaluate a local triviahzation : 7r~^{U) — > [/ x R'^ of the 
vector bundle V as follows: 

To)) = (x, {X^(x),u^(x, ro))x, {X'^(x),u^(x, ro))^), F e F^,^^. 

Let f/a n f//3 7^ and A^, ■ ■ ■ , X'^ and Xj^, . . . , A^ be frame fields over Ua and f/g, 
respectively. Then for any x e Ua r\Uf}, Xa{x) , . . . , Xa{x) and Xp{x), . . . , A^(x) 
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are bases of tt ^{x). Thus, any vector u^{x,to) G tt ^{x) can be expressed as 

^ d ^ d ^ 

i=l 'i=l 

Therefore, one can find a matrix G{x) e Gh{d, M.) such that 

{{Xl{x),u^{x, To))x, • • • , (A^(a;), M^(a;, to))x)G{x) 
and hence, the GL((i, M)-cocycle tapix) = G{x). 

□ 

Lemma 4.23 and Theorem 2.1 imply that the bundle is reconstructed. 

4.4. Reconstruction of the Riemannian structure on the vector bundle. 

Let Jq, /i e be some intervals such that Ii < Iq and Tq e Iq. Then consider 
the set 

5 = {F e :f : To) = for all x e M}. 

It is clear that S ^ $. 

Lemma 4.24. The set S can be computed using admissible energy measurements. 

Proof. Let 0a : 7r~^(C/a) — x R'^ be an arbitrary local trivialisation of the bun- 
dle V and A^, . . . , A^ be a frame field over Ua which we have already determined. 
Then F e S ii and only if 

{u^{x, To), X''^{x))x — 0, for all x e Ua, I — 1, ... ,d, and all a 
which can be verified using admissible measurements. □ 

Now let X E Ua he arbitrary point and consider B{x,r) C Ua- We use the 
following formula to recover the inner product of the wave dtu^{x, Tq), F e 5, 

As the manifold and measure on it are known, we can find vo\{B{x, r)). 

Now let us explain how to determine the L^-norm of the wave dtU^{To) in the ball 
B{x, r). First of all since F e <S, by definition of energy, we can find L^-norm of 
the wave dtU^{ro) on the whole manifold M, 

\\dtU^iro)\\h(^M) = 2 lim lim F^(i^^ro) = 2Ea{F,to), F = ((F,^)~ ,)r=i- 
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Then let 

K = K{F,x,r) = {H eS: dtu"{x,To) = dtU^{x,To) for all x G B{x,r)}. 

The set K can be determined using energy measurements. Indeed let A^, . . . , Aj^ 
be a frame field over Ua which we have already found. Then H & K ii and only 
if 

{dtu^~^{x, To), \^^{x))x = 0, for all x e B{x, r), I = 1, . . . ,d. 

By Corollary 4.14 the above conditions can be verified using admissible energy 
measurements. 

Next 

Hu^i^o)\\l2(B(x,r),V\s^s,r)) = \\9tu"i^o)\\l2(M,vy 

Thus, using (4.15) we can determine inner product of the wave dtu^{x, tq), F & S, 
with itself, for all x E Ua- 

Now 

7r~^{x) = {X{x) : A is a C°° smooth section} = {dtU^{x,To) : F e S}. 

Hence we can recover the inner product {X{x), X{x))x for any X(x) e t:~^{x). 
Using the polarization formula for a real vector bundle, we can find 

(A(x),A(x))^= ^(||A(x) + A(x)|||- ||A(x) - A(x)lll), 
for any X{x), X{x) e 7r~^{x), x e Ua- 

4.5. Reconstruction of operator A. To reconstruct an operator A on the 
vector bundle V it is enough to find its representation in an arbitrary local triv- 
ialization (pa ■ 7r"^(t/a) Ua x M"^. 

For any F e JT^^ by Lemma 4.15 and Corollary 4.16 we can find the represen- 
tations of u^{x,To) and dfU^{x,TQ) in the local trivialization. Thus, using the 
wave equation we can find 

Au^{x, To) = -dfU^{x, To) 

in the local trivialization. Recall that any v e H^'^^{M,V) can be written as 
V — u^(ro) with some F e •^7i,ro- Thus, the set of all pairs 

{{u^{ro)\u.,-u'?^{x,To)\uJ:Ferj^^^J 

is the graph of the operator A in this local trivialization. This determines the 
operator A on M. 
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